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For a symmetric bridge coupled to infinite leads, in the presence of a dipole-coupled external 
ac-field with harmonic mixing, we solve the Schrodinger equation in the time-domain using open 
boundary conditions as well as in the energy-domain using Floquet scattering theory. As this 
T-H ' potential breaks parity and generalized parity, we find a non-vanishing average current. We then 

T-H \ optimize the relative amplitude ratio between the fundamental and the second harmonic leading to 

c 2 ^ , a maximum in the pump current. 

(N ■ 

'. I. INTRODUCTION 

' Generating non- vanishing average electronic currents in statically unbiased systems, i.e. electron pumping, by exter- 
\ nal ac- fields has been realized in experiments with quantum dots [l|-|3| , nanotubes [4| , semiconductor heterostructures 
[1,0] and a Josephson junction array 0. 
I Theoretically, the adiabatic case of slow driving has been treated amongst others by Thouless by Brouwer 

. 0, who used a scattering approach, and by Zhou et al. [13], applying Keldysh-Green's function methodology. 
(-H ' The non-adiabatic, time-periodic case can be handled in the Floquet formalism, formally treating the system as 
I , time-independent. An early example of a similar strategy is the heuristic approach to understand the effect of time- 
^ ' periodic driving on the current vol tag e characteristics of superconductor-insulator-superconductor junctions that has 

a been given by Tien and Gordon lllll . More recently, Floquet scattering theory has been used to get expressions 
for transmission probabilities [l3 - [l8{ . In addition, Kim [l3] as well as Moskalets and Biittiker f2Q| exploited this 
method for considering an electronic pump consisting of two oscillating ^-peaks. In [l^ - fTsl . [T7l - l2Clj ] the scattering 
matrix is determined from the matching conditions of the wave function. Furthermore, electron-pumping scenarios 
have been treated with the help of Floquet theory using an equation-of-motion approach in the Heisenberg picture 
l^^], solving the master equation '22U24'], or using non-equilibrium Green's functions [25l - l27| . The authors of 
2^,.26j have considered tight-binding Hamiltonians, while in pSj an interacting two-level system has been studied. 
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21, 
21- 

27j provides a comparison of Green's function theory with the scattering matrix. 

Recently, Kurth and coworkers [285 proposed a time-dependent approach for quantum transport, treating the leads 
by using open boundaries for the central region. This method is not restricted to time-periodic problems. In the 
fS| same Reference a simple scheme for ac-transport is studied as one application. Stefanucci et al. [29j later used this 
^ ' algorithm of time-evolution for the investigation of systems where a net current is generated by a traveling wave in 
, the potential of the Schrodinger equation. 



In this paper, we consider a harmonic-mixing dipole field as another potential that shows the pumping effect. In 
order to break (generalized) parity, the external field is composed of a fundamental frequency component together 
with an additional second harmonic term. Harmonic mixing has been studied previously in the tight-binding case 
<N ■ f2l[2l as well as purely classically [sol]. In Section |IT1 in order to calculate the transport across a structure under 
\^ . the infiuence of the external ac-field, we first review Floquet scattering theory, allowing us to arrive at an expression 

for the stationary net current from the matching conditions via the scattering matrix. 
• • 1 Secondly, in Section IIIIl we treat the transport problem from a time- dependent point of view, focusing on the 
_ ^ ' transient dynamics of the current after a sudden switching on of the driving. To this end we are considering the zero 
K> ', temperature case and are choosing an equilibrium state of the undriven system as initial condition [sij . We review 
j_j ' the algorithm of time-evolution from Ref. [28| , where an expression for the time-dependent current is gained by using 
I a Schrodinger equation with open boundaries. In Section |IV] we explicitly apply both ways to calculate the current 
in the problem of electron pumping and it is shown that the time-average of the current over one period converges 
in the long-time limit to the result of Floquet theory. That was to be expected as there are no bound states in the 
considered system which could lead to a non-convergent time dependence as shown in the work of Khosravi et al. 
[33 ]. Our studies extend previous work in the monochromatic case by Li and Reichl 15] and lead to optimal values 
for the relative amplitudes of the first and second harmonic as a function of the incoming energy. In Section |V] we 
give conclusions and an outlook. 
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II. FLOQUET SCATTERING THEORY 



Floquet scattering theory has been employed to calculate transmission probabilities in tunneling systems driven by 
a monochromatic laser field [isl . [33| . Here we briefly review the formalism in order to apply it to the phenomenon of 
electron pumping. 

We consider the time-dependent Schrodinger equation (TDSE) with one dimension in space (atomic units (a.u.) 
are used throughout the paper), 



-ldl + V{x,t) 



^ix,t). (1) 



The potential shall fulfill the following conditions: 

1. V{x,t + T) = V{x,t) (periodicity in time), 

2. there are xl < xr £M. with V{x, t) ^ Q li x < x^ ot li x > xr, 

3. within xl < x < xr we know the analytical solution of the TDSE, 

4. the potential is bounded. 

We denote the left lead (x < x^) with L, the right one {x > xr) with _R, and the central region (xl < x < xr) with 
C. Due to the time periodicity, i.e. condition 1 above, solutions of the TDSE can be written in the form [g^] 

^{x,t) = e-'''^ix,t), (2) 

OO 

cj){x,t + T) = cl>{x,t)^ J2 c„(a;)e-™'^*, (3) 

n— — oo 

where e denotes the so called Floquet energy or quasi-energy, the frequency oj is given by = 27r/T. 
The second condition for V{x^t) leads to the following form of the solution in the region x < xl-, 

*(x, t) = e-'"' Yl ( -^e*'^"^ + ^e~''="^ ) e"™"* (4) 

and for x > xr, 



with the wavenumbers 



k„ = y^2{e + nu}). (6) 
The Floquet scattering matrix S connects the outgoing current amplitudes b^^ with the ingoing ones a^: 

K= J2 E'^W- (7) 

fl=L,R m 

To quantify the net current (/) = 1/T Jq dt I{t), we have to include a summation over the Floquet modes in the 
Landauer formula [l^l as an incoming wave with energy E can be scattered to an outgoing one with energy E + nuj 
(n G Z). In the zero temperature case and including a factor two for the spin this leads to 

W - ^ / E [\STl{E)\' - \Slt{E)\'] ^ JdE^. (8) 

Q m>mo Q 

Here the sum is only taken over those m which fulfill m > mo with mo = —E/uj. That means bound states, for which 
kn is imaginary, do not contribute to the current. 

\S'^l{E)\'^ = |S'rfl(-E)P holds and thus (/) vanishes if the potential fulfills parity {V{-x,t) = V{x,t)) ^ or 
generalized parity {V{—x,t + T/2) = V{x,t)) [l^. As we are interested in non-vanishing net current the potential 
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to be studied below is supposed to break (generalized) parity. This can be achieved either by a breaking of the 
symmetry in position space (ratchet effect) and with an additional monochromatic field or by temporal symmetry 
breaking (harmonic mixing) [2^. We will consider the second case. In order to calculate the relevant part of the 
Floquet scattering matrix that appears in Eq. we use matching conditions for the solutions and their derivatives 
at the boundaries of region C . More details for the case to be considered in Sec. IIVI can be found in Appendix A. 

The formalism can be generalized for cases in which the analytical solution is not known for the whole region C 
but for sections composing C and for cases in which the potential in the leads is not zero but position- independent. 



III. TIME EVOLUTION 



Discretization of the TDSE in space (interval length Ax) leads to a tridiagonal matrix for the Hamiltonian, 







n 


H = 


n 


hi n 







n hi+i 



(9) 



with 



and 



n ■ 



2{Ax) 



(10) 



Note that H = H{t) is explicit time-dependent 21 By splitting ip into the projections on the three regions x < xl, 
X £ [xLjXfl] and x > xr and calling -0 in the respective region tp^, ipc and ipu the TDSE takes a form which has 
been studied by Heliums and Frensley [si" 




Hll Hlc \ / V^L 
HcL Hcc{t) HcR i>c 

HrC HrR / \ t^R 



(11) 



Haa (q^ — L, C, R) are tridiagonal matrices and Hca, Hac {oi — L, R) have only one non-zero entry. We restrict our 
discussion in the remainder of this paper to the case of an explicit time dependence of the Hamiltonian only in the 
central region. 

The idea of handling (jlip is to treat only the central part x £ [xl, xr] explicitly as an open quantum system and 
get the influence of the leads (that are regions x < x^ and x > xr) by a source-term giving the influence of the wave 
function in the leads and a memory-term giving the feedback of the part of ipc(t = 0) which propagates into the leads. 
This has been done in for a static Hamiltonian by calculating expressions for the propagator in the leads. 

Kurth et al. (28j presented a numerical scheme for solving equations given in the form (jlip even for cases with 
a time- but not position-dependent potential in the leads using a generalized Cayley method. In the case of zero 
potential in the leads, it has the form 



(1 + i5i?('"))V'(™+i) = (1 - i5ij(™))V'^™\ 



(12) 



wherein (m) denotes the index of time, = i (H{t.m+i) + H{t„i)), — mAt and S — At/2 is a half time step in 

the time discretization. Splitting (I12p into parts according to the spatial segmentation and employing the parts for L 
and R in the one for C leads to 



[1 + i5Hlp]4 



a=L,R 



with i?^^'' and Ta™' defined as 



h!,P (i^cc^ - iSHcLil + iSHLLV^HLc - iSHRc{l + tSHRR)-^HRc 



(13) 



(14) 



1 — iSHaa 
1 -I- iSHaa 



(m) 



(15) 
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As V'q"'* will not be calculated, it has to be replaced in dTSl) via the projection of (fT2l) onto L, respectively R. The 
result reads 



2iSH, 



Ca 



(1 - iSHg,^)^-'' 

1 V-^a ■ 



fe=l 

(1 - iSHaa) 



(1 + i5H^„)--+^ 



(1 - 

(1 + i5H^^)^-k+-i 



(16) 



The first line of is called memory term as it describes the influence of the history of tAc up to the time tm- 
The second line gives the effect of the initial wave function in the leads and is called source term. In [28] a detailed 
instruction for implementing the equations above is given. In this work we follow those lines. 

At t — we start in the ground state of the total system without the time-dependent part of the potential, which 
is switched on at < = 0. For the calculation of the current, we have to take into account all values of the wavenumber 
k between zero and kp = \/2Ep where Ep denotes the Fermi energy. For each k there are two linear independent 
wave functions. In the leads they are given by Eq. (|4l[5]) with = 0, = 5on (denoted as ij}^) and with = (5o„, 
= (denoted as V'l')- The ^-coefficients follow from the (standard) scattering matrix for i < 0. The initial wave 
function in the central region are obtained from the matching conditions in a similar way as the scattering matrix. 
We get the current by integration over k [2^ : 



iix,t)^ I dk is ((V'^)*a,^,^ + (V'|^)*9,vf ; 



(17) 



dl/dk 

To compare this to the calculation with the Floquet scattering matrix we use that d//d£' = d//dfc x 1/fc. 
As it is done in we define the temporal mean value of the current (/) {x, t) (for fixed position x) by 



{I){x,t) 



Q{T-t)- f dt' I{x,t') + Q{t-T)l- j dt' I{x,t') 



(18) 



As a typical behavior of an open system without bound states we expect that deviation from the solution of Floquet 
scattering vanishes for t — > oo and thus lmit^oo{I){x,t) = (/) with (/) from ([8]). Furthermore, we stress that the 
time-dependent formalism with open boundaries goes beyond standard wave-packet approaches |35l44ll | , yielding only 
transmission probabilities and not allowing for a description of the transient time-evolution of the composite system 
dynamics. 



IV. HARMONIC-MIXING DIPOLE FIELD 



In the remainder of this paper we study the transport induced by a dipole field in the central region of our system. 
The potential in the leads is assumed to be zero and the generalized parity is broken by harmonic mixing according 
to either 

V{x,t) = e(t)e(d/2- \x\)x{Aain{ujt) + Bcos{2Lut)) (19) 

or 

V{x,t) = e(t)e(d/2- \x\)x {Acos{ujt) + B sm{2ujt)) . (20) 

We refer to the potentials (fT9|) . (|20)) as case I, respectively II, and will concentrate on the optimization of the relative 
strengths A and B of the first and second harmonic later-on. 

The analytical solution of the TDSE within |a;| < d/2 and a scheme for extracting the Floquet scattering matrix S 
using the matching conditions is given in Appendix \X\ for (ITOl) . Case II of Eq. ([20]) can be treated accordingly. 

As can be seen in Figure [T] there is a qualitative difference in the potentials described by (fTO]) and (|20l). Especially 
time-reversal parity {V{—x,—t) = V{x,t)) is present in ([SO)) after a time shift of tt/{2uj) while it is broken in (|19p . 
where time-reversal symmetry {V{x, —t) = V{x, t)) is fulfilled after t ^ t + ir/ (2a;). 
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A. Asymptotic average current 



Although cases I and II have different symmetry properties, our numerical results for the integrand of the current, 
d{I)/dE in ([8|), in case of a ratio of A/B ~ 2/1 are of the same order of magnitude, as can be seen in Figure [2] 
In the weak-coupling limit it has been found for the tight-binding scheme that the current vanishes in linear order 
when time-reversal parity is present .22„ ,23J . Clearly, there is no weak-coupling limit within the TDSE we solved here 
and therefore in our case II, by effectively going to arbitrarily high orders, we have found a non-vanishing current. 
Furthermore, we stress that the results from Floquet scattering theory and the converged time-dependent ones shown 
in Figure [2] coincide within numerical accuracy. 

The characteristic jumps in d{I) /dE{E) at E = nu! {n G N) displayed in Figure [5] arise from the fact that at those 
energies another scattering channel is opening because E — nuj becomes bigger than zero. In other words, there is 
another Floquet mode with real wavenumber. The sums in ([8]) contain one more term, but that is not the only reason 
for the changes in d{I)/dE because then we would find a step-like change. A characteristic behavior near E = nuj in 
the sense of high or discontinuous slopes can already be found for 15*^^ as a function of E (not shown). The reason 
is that the matching conditions are sensible to changes in E when the wavenumber for a Floquet mode is near zero 
and changes from imaginary to real for increasing E. In the insets of Figure [5] the convergence characteristics of the 
time-dependent results near E = oj are displayed and it is shown that increasing the final time improves the results 
towards the Floquet result. Nevertheless, at i? w w, the convergence is comparative slow. 



B. Time-dependent current 



To show the time dependence of the current after switching on the potential in the cases I and II, {I){x,t) is 
plotted on a logarithmic time scale (Figures [3] and HI respectively). As in [2^ the results for the middle {x = 0) 
as well as for the left and right sides {x — zLd/2) of the central region are shown with different colors (line styles). 
Although I{x,t) seems to be nearly time-periodic after only one period of time (see insets of Figures [3] and S]) , there 
are changes in I{x,t) for longer times, which can be seen in the behavior of {I){x,t) that needs about 10'^ a.u. to 
coincide approximately with the result from Floquet scattering theory. 

Furthermore, we note that for both potentials, the amplitude of the current (see insets in Figures [3] and S]) is much 
higher at x = than ai x — ±d/2. This means, in these dipole fields, there exists periodic forward and backward 
electronic transport which does not arrive at the leads. 
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FIG. 2: (Color online) d{I)/dE for case 1 (a, b) and case II (c, d) calculated with Floquet scattering theory (black solid 
lines) and time-dependent solution (red dashed lines). Parameters are in both cases d — 3 a.u., A = 2B — 0.25 a.u., 
LJ = 27r/15 a.u. « 0.419 a.u. Numerical parameters of time evolution are Ax — 0.01 a.u. and 5 — 0.05 a.u. The time-dependent 
result is plotted for t = 5000 a.u. and x = (the green dotted lines in the insets show the result for t — 10000 a.u.). The 
scattering matrix has been computed with 11 Floquet modes for (a), (c) and 17 for the insets. 



C. Optimization of mixing parameter 

So far we had used the ratio A/B = 2/1 for the relative strength of first and second harmonic. In order to find the 
optimal mixing parameter B^/{B^ + A^) we have calculated the average current ([S]) for a fixed Fermi energy of 
Ep — 0.3 a.u. and for a fixed value o£ + B^ — 0.078125 a.u. using Floquet scattering theory. We refrain from doing 
time-dependent calculations here as the agreements of the results has been shown above and as the time-dependent 
calculation is numerically more costly. One finds that to achieve high values of | (/) | it is not optimal to use a ratio 
of A/B = 2/1 (thus B'^/{A'^ + B^) = 0.2) as done in the numerical examples before. The results presented in Fig. [5] 
shows a minimum (i.e. a maximum in the absolute value) in the mixing parameter at 0.34 for both cases. 

In order to elucidate which energy contributes dominantly to the current, we look at the differential expression 
d{I)/dE defined in ((8]) for different values of the incoming energy E and the parameters A and B in p9|l and (f20|) 
again for a fixed value of A^ + B^ = 0.078125 a.u. In the limits of either i? = or A = 0, d{I)/dE is zero as 
generalized parity is valid in these cases. In Figure |6] the results of the calculation are shown in a contour plot. The 
maximum and minimum of d{I)/dE as a function of B'^ /{A'^ + B"^) for fixed E are highlighted. For several values of 
E an extremum of d{I)/dE as a function of B^/{A^ + B"^) lies around 0.33 thus B"^ /A^ « 1/2 especially in case II. 
On the other hand there are no sharp maxima or minima for variable B in contrast to d(/)/d£' as a function of E 
for a fixed potential. Thus |(/)| depends more sensitive on the Fermi energy than on the mixing parameter. 

Finally, it is known from classical [13] as well as from tight-binding calculations [l^l that typical harmonic mixing 
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FIG. 3: (Color online) Time-dependent results for {I){x,t) (a) and I{x,t) (b) in case I a.t x = —d/2 (blue dotted), a; = (red 
solid) and x = d/2 (green dashed). Parameters are the same as in Figure[2l Ef = 0.3 a.u. and we have discretized the integration 
in Eq. (|17|l for 100 values of k. For comparison the scattering matrix has been computed with matrices with 11 Floquet modes 
and for 2000 energy- values between and Ep- Discretized integration according to (|8]) gives (/) = -1.11 X 10"^ a.u. (black 
lines). 



signals for a time-dependent force F cx ^cos(a;i) + B cos{2u!t + (/>) are in lowest order proportional to A^Bcos{cj)). 
Our cases I and II correspond to cj) = tt (I) and (p = — 7r/2 (II), respectively. The proportionality to A^B is equivalent 
to a proportionality to ^/x — Var^ in the mixing parameter x :— B^/{B^ + A^) for constant + B^, which leads 
to an extremum at x = 1/3. For the situation considered here, the numerical results show that this dependence on 
the mixing parameter can be found for a large range of incoming energies (see Fig. [S]) and in the integrated current 
plotted in Fig. [5l The simple dependence on the phase 0, however, is not observed in general. 

V. CONCLUSIONS AND OUTLOOK 

We have investigated a potential with broken (generalized) parity in order to generate non-zero net currents in 
unbiased systems. As expected, in the long-time limit, the values from time-dependent calculations using open 
boundaries agree with the results of the Floquet scattering theory obtained from the matching conditions. The 
characteristic behavior of d{I) /dE{E) where E/uj is an integer can be understood as the opening of another Floquet 
channel for the scattering. Optimization of the mixing parameter _B^/(A^ -t- B^) shows that for a large range of 
energies a value close to 0.33 leads to a maximum current amplitude. 

In our example the potential is zero in the leads but the applied methodology can also be used for systems with 
time-dependent potentials in the leads. Floquet scattering theory is numerically less costly than time evolution for 
the case considered here, but it is restricted to problems where the time dependence is periodic and the analytical 
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FIG. 4: (Color online) Time-dependent results for (I) (a) and / (b) in case 11 at a: = —d/2 (blue dotted), a; = (red solid) and 
X — d/2 (green dashed). All calculations has been done for the same parameters and in the same way as for Fig. [S] The result 
of Floquet scattering theory is (J) = —1.21 x 10^"^ a.u. (black lines). 



solution of the TDSE must be known within certain intervals. Moreover a numerical treatment is only possible if a 
finite number of Floquet modes have a substantial contribution to the full solutions of the TDSE. 

The time-dependent calculation, however, is not restricted to periodic problems and it provides information about 
the time it takes after switching on a potential to arrive at the quasi-stationary Floquet results. This time may be 
of interest for microelectronics applications. The algorithm of time evolution can be used for tight-binding problems 
as well as for the continuous TDSE considered here. For the solution of Kohn-Sham equations within TDDFT the 
algorithm has been used recently by Kurth et al. 42]. 
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Appendix A: Floquet S matrix for the harmonic-mixing problem 

In this appendix we give details of the Floquet approach to scattering in time-periodic potentials for the harmonic 
mixing problem studied herein. 

To start the discussion, we first briefly review the case of a monochromatic potential of the form V{x, t) = Ax cos{ujt). 
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FIG. 5: (/) as a function of B'^I{B'^+J^) for case I (solid) and case II (dashed). The energy integral in ([U has been discretized 
by using 100 energy values between and Ef = 0.3 a.u. Calculations have been performed again with 11 Floquet modes. 



There the so-called Volkov-solution of the TDSE reads [15[ 

oo 



n— — oo 



wherein the definitions 



/■* A 

f(x,t) = et+ dt' V(x,t') + h(t) ^ et + —xsm(ujt) 5-sin(2a;i), 

J w 

g{x,t) = x+ dt' / dt" V{x,t")/x = X ^cos(wi), 



= V2(e- A2/(4cj2)+na;) 



(Al) 



(A2) 

(A3) 
(A4) 



have been used. 

The Ansatz (IA1[) is now made for the harmonic-mixing potential V{x,t) = x [Asm{ujt) + _B cos(2a;t)], with the new 
definitions 



f{x,t) ^ et+ dt'V{x,t') + h{t) = et + : 



cos(wt) H sm(2ajt) 

u) 2uj 



hit), 



g{x,t) = x + f dt' J' dt"^ 



„V{x,t") A 

— X — 



B 



siniujt) ^ cos(2a;t) 

^ 4w^ 



(A5) 
(A6) 



and Qn — y^2(e + Co + nw) with constant Cq. For h{t) after some algebra, we get the following ordinary differential 
equation, 



h = —— cos{2Lut) - ——- cos(4a;t) - — — [sin(3a;t) + sin(a;f)] + — — + — — + Co. 



(A7) 



The choice Cq = —A^/{Auj'^) — 5^/(16^^) warrants that h{t) is periodic so that e is the Floquet energy. Simple 
integration leads to 



A^ AB AB 

^^'^^ " 8^ '"'^^^^^ " 64^ "'"'^^"'^^ + 12Z;^ cos(3c^i) + ^ cos{ut), 



(AS) 
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FIG. 6: (Color online) Contour plot of d{I)/dE as a function of the incoming energy E and the mixing parameter / {A^ + B^) 
with fixed + B'^ = 0.078125 a.u. for case I (a) and II (b). The maximum (black dashed) and the minimum (white solid 
lines) of d(J) /dE for fixed E are highlighted. We take once more Floquet modes with energy between E — 5uj and E + into 
account. 



which completes the solution in the central region C in case I of Sec. IIVI Regarding the potential (|19p for t > 0, 
the analytical solution of the TDSE can be written down in each of the regions x < —d/2 (L), \x\ < d/2 (C) and 
X > d/2 (R), 



n 

Mx,t) = e-"*^ 

n 



L „ik„x 



„ig„g(x.t) 



'iq„g(x.t) 



(A9) 
(AlO) 
(All) 



/9n V?" 

with q„, f{x,t), g{x,t) as determined before and fc„ = •\/2(e + nuj). The matching conditions at xl — —d/2 and 
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xb. = Ajl are 



i>L{xL-^A = ■ipc{xL + 0,t), 

V'fl (a^fl + 0, = Vc (a^fl - 0,i) 7 
V'i? (a^fl + 0, i) ^ ■il;'c{xR-0,t). 



(A12) 
(A13) 
(A14) 
(A15) 



Employing the operator ^ Jq"^ dt e^^^^ {s G Z) on the matching conditions above, i.e., make a Fourier transformation, 
results in 



n ijn 



sn 



(A16) 



iqnd/2 



71 ^ 



(A17) 



E 



-iq„d/2fjR- 



(A18) 



-iksd/2 



x<i/2 



_|_^"g-i9nd/2 



sn 



where the definitions 



a± 



l/T I dt exp(i [(s - - ± 9nG(t)]), 



Z?^± = / di F'(i)exp(i[(s-7i)cjf-F(x„,i)±9„G(<)]) 

Jo 



(A19) 



(A20) 



(A21) 



have been introduced with a — L,R and F{x,t) ~ fix,t) — et, G{t) = g{x,t) — x, T = 2tt/uj. By introducing the 
diagonal matrices {E)km = (5fcm exp(?gmd/2), {X)krn = (5fem exp(ifcmd/2), (Q)fcm = (5femV9^) (-?^)fcm = Skm\/%^, it is 
possible to write the matching conditions in matrix form, 

K-^X-'^UL+K-'^XbL = G^+Q-^E-^a + C^-Q-^Eb, (A22) 

KX-^QL-KXbL = [C^+QE-^ +D^+Q-^E-^]a 

+ [-C^-QE + D^-Q-^E] b, (A23) 

C"+Q-^Ea + C^-Q-^E-^b, (A24) 

[C"+QE + D^+Q-^E] a 

+ [-C"-QE-^ + D^-Q-^E-^] b. (A25) 

To get an expression for the Floquet scattering matrix S defined in Eq. ([T]) we need b^ and b^ in dependence of 
gl and ur. By multiplying Eq. (jA22l) from the left with the matrix K and Eq. (|A23p with it is possible to 
eliminate the coefficient b^. For we find the following expression in dependence of a and b, 



R-^X-'^aR - KXbR 
-KX-^GR + KXbR 



■Bi 



y [RC^-Q-^ - R-^C^-Q + R-^] E b. 



(A26) 
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In the same manner we can eliminate bn from Eq. (jA24l) and Eq. (jA25|) and get 



^ (A27) 

+ — [KC^-Q-^ + R-^C^-Q - R-^D^-Q-^] E'^ b. 
. ' 

Eliminating now the coefficients for the incoming waves and a/j from Eqs. (jA22|) . (jA23|) or (jA24[) . (jA25|) . respec- 
tively, we find 



bL=^ [rC'+Q-' - R-'C^+Q - R-'D'^+Q-'] a 

" V ' 

= :B5 

+ — - [RC^-Q-^ + R-^C^~Q - R~^D^-Q-^] Eb 



(A28) 



and 



bn^^ [RC"+Q-' + R-'C'^+Q + R-'d'^+Q-'] E a 



= :B 



(A29) 



[RC^-Q-^ - R-^C^-Q + R-^D"-Q-^] E-^ b. 



=-Ba 



Converting Eq. (|A26|) leads to a = ^ (a^ - B2b). Putting this in Eq. (|A27p results in BsB^ ^ (ol - -82^) + B^b = 
qr. Converting that to an expression for b yields 

b=[Bi- 333^^82] {an - BsB^'aL) . (A30) 

One the other hand one gets from Eq. (|A27P b = B^^ (ajj — B^a). Putting this in Eq. (jA26P then leads to 
(Bi — 52-847^-83) a + B2B^^aB. = cll and thus wc find 

a=[Bi- S2S4 (flL - -82^4-10^) . (A31) 
V ' 

Putting the equations (|A30|) and (|A31I) in Eqs. (|A28|) and (IA29P we directly get 

bL = SsAfi-i (ttL - B2B^^aR) + BeM2^ (an - BsB^^aL) , (A32) 

6i^ = ByM-^ (fli - B2B-^aR) + BgM-^ {ur, ~ BsB-^ul) . (A33) 

We find 

Slr = -BsMf IB2-84-1 + S6A/2-\ (A34) 

Srl = BrM^^-BsM^'BsB^K (A35) 

As only Slr and Srl are needed to calculate (/) in Eq. we have found expressions for the relevant parts of the 
Floquet S matrix for the considered system with harmonic mixing. 
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